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Abstract. The diffusive thermal conductivity of neutrons in dense matter [p ~ (1 — 8) x 10^* g cm"''] of 
neutron star cores is calculated. The contribution from neutron-neutron and neutron-proton collisions is taken 
into account. We use the transition probabilities calculated for symmetric dense nucleon matter on the basis of 
the Dirac-Brueckner approach to the in-medium effects and the Bonn model of bare nucleon-nucleon interaction. 
The diffusive thermal conductivity of neutrons in the presence of neutron and proton superfluidities is analyzed 
in a microscopic manner; the effects of superfluidity are shown to be significant. The low temperature behavior 
of the thermal conductivity appears to be extremely sensitive to the relation between critical temperatures of 
neutrons and protons. The results are fitted by simple analytic expressions. In combination with the formulae 
for the electron and muon thermal conductivities, obtained earlier, the present expressions provide a realistic 
description of the full diffusive thermal conductivity in the neutron star cores for normal and various superfiuid 
phases. 
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1. Introduction 

We study the thermal conductivity of neutrons in the neu- 
tron star cores, in the density range p ~ (0.5 — 3)poi where 
po = 2.8 X 10^^ g cm^'^ is the normal nuclear density. The 
problem is important for numerical simulations of cooling 
of young neutron stars, before internal thermal equilib- 
rium is achieved. The non-relaxed stage lasts 10 — 100 
yr (Van Riper, 1991, Umeda et al., 1993, Lattimer et al., 



1994) depending on the model of dense matter. 

At densities close to po, neutron star matter is ex- 
pected to consist mostly of neutrons (n), with a small (a 
few percent) admixture of protons (p), electrons (e) and, 
at densities where the electron Fermi energy exceeds the 
muon rest energy, also of muons (/i). In what follows, we 
adopt this npep, model of neutron star matter within the 
considered density range. At higher p other particles may 
appear, first of all hyperons (e.g., Balberg et al. 1999). We 
discuss briefly thermal conductivity of hyperonic matter 
in Sect. 4. 

Numerical simulations of neutron star cooling tra- 
ditionally employ the thermal conductivity obtained by 
Flowers & Itoh ( |1979| , |198l| ) about 20 years ago. These 
authors proposed an analytic fit valid for non-superfiuid 
matter, and discussed qualitatively the effects of nucleon 



superfluidity. It is now generally recognized that the neu- 
tron star cores can be in a superfluid state (see, e.g.. Pines, 

Yakovlev et al. 



1991, Page & Applegate, 1992 



Moreover, Flowers & Itoh ( 1979| , 



1999) 



1981) calculated the ther- 
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mal conductivity for one specific model of dense matter 
which is inconvenient as the equation of state of neutron 
star cores is generally unknown. The thermal conductivity 
of the neutro n star cores was analyzed a lso by Wambach 
et al. (1993) and by Sedrakian et al. ( |1994|) , again for 
some selected models of dense matter and neglecting the 
superfluidity effects. Thus, it is desirable to obtain new 
expressions which apply to various models of matter and 
take into account nucleon superfluidity. In contrast to the 
studies of Flower s fc I toh ( |l979| ), Wambach et al. ( [1993D , 
Sedrakian et al. (1994), we consider the effects of nucleon 
superfluidity in a microscopic manner by evaluating the 
diffusive thermal conductivity of neutrons in superfluid 
matter. 

The main heat carriers in the neutron star cores are 
neutrons (the most abundant particles) as well as electrons 
and muons which have large mean free paths since they ex- 
perience only relatively weak. Coulomb interactions. The 
therm al conductivity of protons is small (Flowers & Itoh, 
1979 ) and can be neglected. The heat conduction of elec- 
trons and muons is limited by Coulomb collisions with e, 
p and /i, while that of neutrons is limited by their col- 
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lisions with n and p due to the strong nucleon-nucleon 
interaction. Accordingly, the neutron transport is almost 
independent of the electron and muon one, and can be 
studied separately. This property was proved by Flowers 
& Itoh (1979) who decoupled the overall heat-conduction 
equations into two blocks describing the conductivities of 
neutrons and electrons. 

The thermal conductivity of electrons and muons was 



reconsidered by Gnedin & Yakovlev (1995) who obtained 
simple expressions valid for a wide class of models of su- 
perfluid and non-superfluid matter. In the present article, 
we reanalyze the thermal conduction of neutrons utilizing 
some new developments in the nucleon-nucleon interac- 
tion theor y. In c ombination with the results of Gnedin & 
Yakovlev ( 1995 ), this describes the diffusive heat trans- 
port in not too dense neutron star cores. 

In the presence of neutron superfluidity there may be 
another channel of heat transport, the so-called convec- 
tive counterflow of normal component of matter with re- 
spect to superfluid one. This mechanism is known to be ex- 
tremely effective in superfluid helium (e.g., Tilley & Tilley, 



tides can be described including only binary quasiparticle 
collisions. We are interested in the collisions of neutron 
quasiparticles with neutron and proton quasiparticles. In 
the calculation of the neutron thermal conductivity pro- 
tons will be considered as scatterers. 

2.1. Kinetic equation 

The distribution function Fn of neutron quasiparticles sat- 
isfies the Landau kinetic equation (e.g., Baym & Pethick, 



1991 ) which describes motion of a given quasiparticle in a 



1960) . The thermal flux 



self-consistent field of its neighbors. However, a proper lin- 
earization reduces the Landau equation to a Boltzmann- 
like equation. The linearized expression for the quasiparti- 
cle thermal fiux density also takes its standard form. This 
enables us to use the expressions familiar from the kinet- 
ics of dilute gases (e.g., Ziman, 
density of neutrons is 

Qn = J dkniCn- Hn)VnFn= -K^S/T, (1) 

^) but in the case of neutron star matter the situation ^^^^^ neutron wave-vector, veloc- 

ity, energy, and chemical potential, respectively; k„ is the 
neutron thermal conductivity to be determined, and T is 
the temperature. In the presence of a weak temperature 
gradient, the Boltzmann equation reads 

where /„„ and /„p are the nn and np collision integrals: 

1/3 



is more complicated. The related effects require separate 
study and will not be considered here. 

2. General Formalism 

Consider a neutron star core at densities p from about 
0.5po to 3po- The lower limit corresponds to the core- 
crust transition (Lorenz et al., 1993, Pethick & Ravenhall, 
1995| ), while the upper limit is the typical central den- 
sity of a not too massive star. We treat neutrons as non- 
relativistic particles and do not study higher p where they 
become mildly relativistic. 

To calculate the neutron conductivity, it is sufficient to 
consider the nucleon component of matter. In the density 
range of study, typical Fermi energy of neutrons changes 
from a few ten to a few hundred MeV (e.g., Shapiro & 
Teukolsky, 198E). Fermi energy of protons varies from sev- 
eral MeV to a few tens of MeV. Since the temperature in 
the cores of cooling neutron stars is typically well below 
10^° K (1 MeV), both neutrons and protons form strongly 
degenerate Fermi liquids. 

In the limit of strong degeneracy transport properties 
of a nucleon liquid can be described using the concept of 
quasiparticles. In the present section, we restrict ourselves 
to the case of normal nucleon liquids (the effects of nu- 
cleon superfluidity will be considered in Sect. ||). Landau 
theory of normal Fermi liquids enables one to identify the 
low temperature properties of a strongly interacting real 
system with those of a dilute gas of weakly interacting el- 
ementary excitations — the nucleon quasiparticles (a de- 
tailed description of the the ory a nd its applications can be 
found in Baym & Pethick, 1991). This formalism can be 
used effectively only for describing small macroscopic de- 
viations of the real system from its ground state. A crucial 
feature of the Landau Fermi-liquid theory is the inclusion 
of the quasiparticle interaction. The kinetics of quasipar- 



j j j Ak[Ak2Ak'^ 



(27r)9(l-t-(5i2) 
X [F[F!,(l-Fi){l~F2) 

- Fi F2 (1 - F[) (1 - F^)] . 



(3) 



The non-primed quantities refer to the particles before a 
collision event, while the primed ones refer to the particles 
after the collision, a denotes the spin state, and 7^12 is the 
differential transition rate. The rates of direct and inverse 
processes are equal due to detailed balancing. The term 
5i2 serves to avoid double counting of the final states in 
the case of nn collisions; V is the normalization volume. 

Since the protons are treated as scatterers (but not the 
heat carriers) their distribution function coincides with 
the equilibrium Fermi-Dirac distribution /(e): 



p — /p — f{^p)j 



1 + exp 



T 



The neutron distribution can be written as 
p — f — ffi 

den 



(4) 



(5) 



where /„ = /(e„), and = <&„(£„) describes the small 
deviation of the distribution function from the equilibrium 
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one. This correction is known (Ziman, 196C) to have the 
form 



X dko dfco 



Vl2 fl /2 



vn ■ vr 

T ' 



(6) 



where Tn is the effective relaxation time of neutrons which 
is generally a complicated function of the energy param- 
eter (e„ — fj,n)/T. However, while solving the transport 
equation, we will rely on the variational approach and 
choose the simplest trial function consistent with a sym- 
metry requirement imposed on <J>„. This corresponds to 
r„ independent of energy. 

Notice that the thermal conductivity can be calcu- 
lated exa ctly u sing the formalism devel oped by Sykes & 
Broo ker (|1970D , H0jgaard Jensen e t al. (|1968| ), Flowers & 
Itoh ( |1979D , and Anderson et al. ( |1987D . In Sect. |^ we 
discuss the corrections to our results which follow from 
the exact theory. 

2.2. Relaxation time 

Inserting F„ from Eqs. and @ into (|^) we have 



3m* 



(7) 



where n„ and m* are the neutron number density and 
effective mass, respectively. The relaxation time r„ has to 
be determined from Eq. (||). 

Following the standard procedure we linearize the ki- 
netic equation with respect to the small nonequilibrium 
correction $„ to the Fermi-Dirac distribution. The lin- 
earized left-hand side of Eq. (||) is 



Vn ■ VFn « Vn ■ V/„ = ^ t;„ • VT. 



(8) 



Inserting Eqs. (jj) and (jg) into (g), we arrive at the lin- 
earized collision integrals 



Il2 = 



(27r)9(l + <5i2)T2 

X (1 - /{) (1 - /^) J2 F 



dfc'i dfc2 dk'2 fl /2 



(9) 



where y„„ = (ei - /i„) t>i + (£2 - M«) '"2 - (e'l - ^J■n) K - 
(4 - Mn) and Vnp = (ei - /i„) Vi - {e[ - v[. 

Let us multiply the linearized Boltzmann equation by 
dkn (27r)"'^ (e„ — Vn, integrate over dfe„ and sum over 
the neutron spin states an ■ From the left-hand side of the 
Boltzmann equation we obtain 



\ 2 ^ .m 



dfn TT^TUn 



(10) 



On the right-hand side, we have 7r2T„T?T,„(^'„„-ft'„p)/m*, 
where we have introduced the effective nn and np collision 
frequencies Vnn and Vnp- 



X {l-f[){^-Qiel-^in)Vl2■Vl. (11) 

Now the relaxation time is: 

rn = ^ . (12) 

The differential transition rate summed over initial and 
final spin states can be written as: 

X S{ki + k2 - k[ - k'^) 



(13) 



where Q12 = (1/4) lisping | A/12 P is the squared matrix 
element averaged over initial and summed over final spin 
states. Since no specific spin orientations are involved, Q12 
can depend only on the scalar combinations of the momen- 
tum transfers q = fej, — ^2 and q' — k'2 ~ ki, that is on 
q, q' and {q ■ q'). In our case, the scalar product [q ■ q') is 
zero (see Sect. ^). Therefore, Q12 is a function of q and 
q' alone. 



2.3. Energy-angular decomposition 

Since neutrons and protons are strongly degenerate, we 
can use the conventional energy-angular decomposition 
while evaluating the collision frequencies V12. The main 
contribution to the integrals comes from narrow vicinities 
of k near k = k-p m which 1iv-p\k — k-p\ ^ T (in this case 
kp and vp are a particle Fermi wavenumber and velocity). 
In these regions, we can introduce convenient momentum 
variable ^ = hv-p{k — kp) with d^ = de. Then we can 
replace integration over k by that over ^ extending the in- 
tegration limits to ±00, since the domains of integration 
outside the indicated vicinities give exponentially small 
contributions. Thus h^dk = kpm* d^ dQ, where df2 is a 
solid angle element in the direction of k. Moreover, we 
can place the particle wave vectors k at the appropriate 
Fermi surfaces in all smooth functions of the momenta. 
Consequently, the expressions for 1/12 are decomposed into 
sums of products containing decoupled integrations over 
angles dCl and momentum (energy) variables d^ of all col- 
liding particles. Placing all the particle wave vectors on 
their Fermi surfaces, one can easily verify that the mo- 
mentum transfers q and q' are perpendicular, as already 
mentioned in Sect. 2.2. Using the well-known symmetry 
properties of the energy integrals and the relations which 
come from momentum conservation for the particles at 
their Fermi surfaces in the angular integrals, after stan- 
dard transformations (e.g., Ziman, |1960D we obtain 

3m^A:FiTO2^fcp2 



1^12 



267r8 T3 {1 + 612) 



da d^[ d6dc^/i/2 



X {l-f[){l-f-)6{ei+e2 



4) 
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dQi dn[ dfl2 dn'^S{ki +k2- k[ - k'2) 

X Qi2(g,g')(ei-ae;^i2), (w) 

where A„„ = 3 - {q^ + q'^)/kl^„ A„p 1 - gV(2fc|,,J. 
These expressions contain two energy integrals Ji and J2 
of the form 

Jc = J J J y^"deideid6de^5(ei+e2-e'i-e2) 

X /1/2 (1-/0(1 -/2)(ei-Mn) Sc., (15) 

with _Bi = £1 — 111 and B2 = — /ii- They may be evalu- 
ated analytically: 



Ji = 3J2 = 

5 



(16) 



Thus we arrive at the equations which contain the angular 
integration 



167r5 mlmf k^^ 

5fcFl(2^)9n^(l + (5i2) 

X <5(fci +k2- k[ - k'2) Qi2(g, q') C12, (17) 



^ iTT m^m^ fcp, r [ f f [dn^ dn[ dn2 dn' 



where Cm ^ q^ + q'^ and C„p = 2fc|„ + (gV2). 
2.4. Angular integration 

It is clear that the integrands depend only on the relative 
orientations of the particle momenta. Then we can imme- 
diately integrate over the orientations of the momentum 
of the first particle before the scattering, and over the 
azimuthal angle which specifies the momentum of the sec- 
ond colliding particle with respect to the first one. This 
yields an overall factor Stt^. Let 9 be an angle between 
the momenta of colliding particles before the scattering, 
and 01, 62, (/>i, </>2 be the polar and azimuthal angles of 
the particles after the scattering in a coordinate frame 
with z axis along ki and x axis placed in a (fci fe2)-plane. 
Accordingly, 



■ a AO ^'^1 
smf^i df^i — -73—, 



sin 62 d92 = 



q' dq' 



(18) 



xpi 



and the integration over 9, 
of the 5-function: 



fcFlfcF2 

'2 is performed with aid 



feF2 

k-pi 



sin 9 d9 ( 



. (5(fci fe2 - fc'i - k'2) 



qq'y/A^-q^ 

where 

= [4fcFl^F2 — (^Fl 



«,p2 



(19) 



(20) 



q< A, l/cFi - fcF2| <q'< (kpi + kF2)- 

Inserting (^9|) into (^7|) we can express Vnn and Vnp in 
a unified manner 



^12 — ; — TT^ '5^i2- 



(21) 



Here, we introduce the convenient quantities of dimension 
of a cross section (cm^): Snn = Sn2 and Snp — Spi + Sp2'- 



Snl 
Sn2 



J I dx' I dx 

Wn'^h Jo Jo 



'AT 



1 



T I dx' / dx 

16n^n^ Jo Jo 



\J\ — — x' 



x^ + 
V 1 — — 



Qnni 



pi 



'AT 



5. 



p2 



167r2a* 



'"■jV 

167r2?i4 



0.5+a;o 



dx' 



a{x') 



dx ■ 



Qr. 



0.5-l-£!;o 



dx' 



0.5— a:n 



a{x') 



dx ■ 



Va2 



Qnp. (22) 



In this case, x — q/{2k-pn) 



= gV(2fcFn), xo = 



(0.25-t-a;^ -a;'2)2/a;', is 



^Fp/(2fcFn), a{x') = 
a bare nucleon mass (we neglect the small difference be- 
tween the bare neutron and proton masses). The quantity 
Snl is not needed here but will be required in Sect. ^. 

Thus we have expressed V12 through the 2D integral 
of the squared matrix element Q12 over the particle mo- 
menta transfers q and q' . This integration can be per- 
formed numerically. Note that for nn collisions q varies 
from to 2A;Fn and scattering at arbitrary angles con- 
tributes to the integration. By contrast, in the np case 
a{x') does not exceed 2xo that is scattering angles larger 
than 2 arcsin(fcFp/fcFn) are forbidden. We recall that, typ- 
ically, /cFp ^ ^Fn, in neutron star cores. 

At this stage it is convenient to introduce the scatter- 
ing cross section. 



2.5. Scattering cross section 

The squared matrix element Q12 determines the scattering 
cross section. The nucleon-nucleon cross section is more 
convenient for practical calculations since it has been the 
subject of some studies. To translate the cross section into 
the squared matrix element one should be careful about 
kinematics. We are looking for a matrix element between 
the quasiparticle states with certain momenta. Let us con- 
sider an idealized case, in which we can approximate nu- 
cleon quasiparticles by bare nucleons. The simplest choice 
would be then to equate this matrix element to that for 
scattering of nucleons with the same momenta in vacuum. 
This proceeds as follows. Let us specify, for our imaginary 
bare nucleons, the center-of-mass (cm) reference frame in 
which the total momentum of a colliding pair is zero, and 
the laboratory (lab) reference frame in which the initial 
momentum of the second particle is zero. Using standard 
quantum mechanics we find the familiar relations: 



dcri2 

dr^cm 



(eiab, 6*0111) 



COS^r 



eiab 



'AT 



167r2?i4 



il2j 



q' 



12 



2mN 



(23) 
(24) 
(25) 
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In this representation, which is popular in the current lit- 
erature, the cross section is a function of scattering angle 
in the center-of-mass frame and collision energy in the 
laboratory frame. 

However, there exists a more involved approximation 
based on recent microscopic calculations of AW scattering 
cross sections in medium (see also Sect. 2.6). To restore 
the squared matrix element in this case, we again con- 
sider the two reference frames but this time for quasipar- 
ticles. The quasiparticle energy and momentum transfor- 
mations which accompany a Galilean transformation from 
one frame to another may be found in Sect. 1.1.3 of Baym 
& Pethick (1991). Although the scattering angle in the 



center-of-mass reference frame is still given by Eq. (p4|), 
other expressions change: 

dai2 ^ (m*m*)2(g2 + g'2)3/2 Q^^ 

X [{k^^ml - kl2ml) {ml - m*) 

X [{q^ +q'^){m\+ml) 

+ {kl^-kl2){ml-ml)]-\ (26) 



eiab 



2mN 



Pf2 ( "^A' 



2m AT \ 7712 



(27) 



If the effective masses of colliding quasiparticles coincide, 
Eq. ( |2^ ) reduces to a much simpler Eq. (p3|), with 77i* in 
place of 777JV. 

2.6. Calculation of Sa 

For the collision energies up to 550 MeV we are interested 
in, the quark structure of matter is not yet pronounced, 
and the bare NN interaction can be described in terms of 
various virtual-meson-exchange processes. 

Very good agreement with experimental data on the 
elastic AW scattering at energies up to 350 MeV is 
achieved using the Bonn potential model (Machleidt et 
al. 1987). It includes all possible diagrams which describe 
w, 77, and S mesons as well as the 
- p mesons and other combinations 



the exchange of tt, p, 
exchange of 27r, 37r, tt - 
of mesons with total energy of any process up to 1 GeV. 
The maximum collision energy of 350 MeV is mentioned 
because the model does not include the inelastic channels 
of TT meson creation allowed at higher energies. We do not 
consider matter with pion condensates, and, therefore, we 
do not allow for these inelastic channels. 

The bare WW potential is known to be attractive at 
large distances and repulsive at short ones. The short- 
range repulsive core is extremely strong, and the Born ap- 
proximation fails to describe the WW collisions of interest. 
Thus we use the exact transition probabilities obtained 
with the Bonn potential. 



The differential cross sections for free nucleons have 
been compiled from tables presented by Li & Machleidt 
These authors consider in-medium effects 



1993 



1994) 



( 

but their calculations are extended to the limit of zero 
medium density which is the in-vacuum case. The latter 
cross sections are computed for laboratory energies up to 
300 MeV, and the results are in nice agreement with ex- 
perimental data. For convenience, we prefer to use these 
theoretical data rather than the experimental ones. The 
cross sections of the nn collisions at higher energies have 
been set equal to their values at 300 MeV since, according 
to the experimental data, the 7771 cross sections are almost 
energy-independent in this range. 

The treatment of the many-body effects is compli- 
cated. For instance, Li & Machleidt (|1993|, |1994[) per- 



formed calculations for symmetric nuclear matter as they 
were interested in heavy-ion collisions. The detailed cal- 
culations for arbitrary asymmetry of matter, in particular, 
for strongly asymmetric matter in neutron star cores have 
not yet been done, to our knowledge. The many-body ef- 
fects turn out to be rather significant; they mainly reduce 
the scattering cross sections by a factor of 2-6. Although 
the in-medium effects can obviously be different in the 
symmetric and asymmetric matter, we make use of the 
available results to illustrate possible influence of medium 
on the transport properties. 

We have calculated the integrals ( p^ ) using both in- 
medium and in-vacuum sets of the scattering cross sec- 
tions obtained with the Bonn model for the values of fcpn 
from 1.1 to 2.6 fm^^ and for the values of fcpp from 0.3 to 
1.2 fm^^ with the maximum ratio of /cpp to /cpn equal to 
0.7. This grid certainly covers the range of possible vari- 
ation of /cpn and fcpp for any realistic equation of state of 
matter at densities O.Spo ^ p ^ 3po- While restoring the 
squared matrix element from the in-medium cross sections 



we set m*^^ m'^^ 0.8 m^ in Eqs. (g6D and (|2 

The numerical results are fitted by simple analytic 
functions and written in the form 



(28) 



where Sa'^ corresponds to scattering of bare particles, and 
Ka describes the in-medium effects. In all these fits, /cpn 
and fcpp are expressed in fm~^. 

The fit expressions for the in-vacuum integrals are 



„(o) 14.57 1 - 0.0788 fcFn + 0.0883 fc2 

S„-{ = —TT- X ; — nrb, 

"1 1- 0.1114 fcpn 

„(o) 7.880 1 - 0.2241 fcFn + 0.2006 H 

Sn-i = -Tt— X ^-nTTTTT — mb. 



1 - 0.1742 fepn 

^(0) ^ 0.8007 fcpp ^ 3^^^g ^ 0.0004285 fc| 



Fn 



Fp 



26.85 fcpn + 0.08012 fc|„) (1 - 0.5898 /cpn 
0.2368 fcpn -I- 0.5838 kl^ -f 0.884fcFn/cpp)"^ mb. 



(0) _ 0.3830 fc^p 

- L.5.5 



(1-f^ 102.0 A:pp + 53.91 fcpn) 



X (1 - 0.7087 fcp„ + 0.2537 k^^ 
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+ 9.404 1.589 fcpn/cFp)"^ mb. (29) 

The mean fit errors of S^^i and 5*^2^ are about 6nns ~ 
0.3%, and the maximum errors (5max ~ 0.5% take place at 
fcpn ~ 1.5 fm~^. In the case of S^p^ , we have Jrms ~ 1.2%, 
with i5max ~ 2.9% at /cpn ~ 1-1 fm""'^ and fcpp ~ 0.3 fm~^, 
while in the case of S^^^ , (Srms ~ 3.4%, with Jmax ~ 9.0% 
at k-pn ~ 1-1 fm^^ and fcpp ~ 0-4 fm^^. 

The fits to the in-medium correction factors Ka are 

Kni = (—\ (0.4583 + 0.8921*2-0.5497^3 

- 0.06205 kpp + 0.04022 + 0.2122 ukpp), 



Kn2 = ( — I (0.4891 + 1.111 u^- 0.2283 

+ 0.01589 /cFp - 0.02099 kl^ + 0.2773 ufcpp), 

Kpi = (^^^ (0.04377+ 1.100^2 + 0.1180^3 

+ 0.1626 fcpp + 0.3871 ukpp - 0.2990u^), 

Kp2 = (^^^ (0.0001313+ 1.2481*2 + 0.2403^3 

+ 0.3257 fcpp + 0.5536 ufcpp - 0.3237u^ 

+ 0.09786 u2/cpp). (30) 

In the case of Kni we define u = kpn — 1.665; the mean 
fit error is about drms ~ 1.3%, and the maximum error 
<5max ~ 3.9% takes place at /cpn ~ 1.1 fm^^ and /cpp w 
0.3 fm~^. For Kn2, we have u = kp^ — 1.556, Srms ~ 1.9%, 
^max ~ 4.7% at fcpn « 1.8 fm""^ and kpp 1.2 fm~\ For 
Kpi, we have u = A:p„ — 2.126, S^-^s ~ 2.4%, Sy^ax ~ 7.1% 
at fcpn « 1.6 fm~^ and fcpp « 1.1 fm^^. Finally, for Kp2, 
we have u = kpn — 2.116, (5rms ~ 3.9%, (Smax ~ 8.5% at 
fcpn ~ 2.1 fm~^ and fcpp w 0.5 fm~^. 

Unfortunately, the quantities Ka depend also implic- 
itly on TO*, as the latter enters Eq. ( p7| ) for eiab. For this 
reason, Eqs. ( pO| ) are strictly valid only if m* = to* = 
0.8TOAr. We will use these quantities for illustration and 
present the final results in the form which could be modi- 
fied easily if more advanced in-medium data are available. 



2. 7. Comparison with exact solution 



As pointed out in Sect. 2.1, we have used the variational 



method while solving the kinetic equation. In principle, it 
is possible to derive an exact thermal conductivity, mak- 
ing use of the advanced theory developed by Sykes & 
Brooker ( 197C| ), and H0jgaard Jensen et al. ( 1968| ), for 
one-com ponen t Fermi-systems, and ext ended by Flowers 
& Itoh (|1979|) and Anderson et al. ( |l987|) to multi- 
component systems. In our case neutrons are the only heat 
carriers, and the equation which determines the energy de- 
pendence of the nonequilibrium term in the neutron dis- 
tribution can be reduced to that for a one-component sys- 



tem: 

+ Xk dx2 



f{x)[l-f{x)]^{x) 



^^±^/(.)/(x2)*(x2). (31) 



In this case, x = {cn — fJ'n)/T, f{x) = l/(e'^ + l), and '^[x) 
is an unknown function defined in the same fashion as in 



Baym & Pethick (1991). The param eter , originally in- 
troduced by Sykes & Brooker ( 1970 ), contains information 
on nn and np scattering and reads 



A 



K 



dn (1 + 2cos( 



1*2 



-Qr. 



COS {e/2) 

j df7[(l + 2x0 cos6')2 +4x0 sin2 6' 



Xo 



C0S( 



(1 + 4x0 cos6' + 4xg)3/2 

^*2 



■xo 



qWLQ 



d» QnnjOA 

cos (61/2) 

n -1 



(1 + 4x0 cos 61 + 4x2)1/2 



(32) 



where 9 and (j) s-rs Abrikosov-Khalatnikov angles, xo = 
fcpp/(2fcFn), and the integrations are carried out over the 
full solid angle (dfi = sin6'd6'd0). In the case of nn col- 
lisions, is the scattering angle in the center-of-mass ref- 
erence frame, while for np collisions we obtain 



(1-4 xg)2 + 16 xg sin^ cos (f> 
(1 + 4xg)2 - 16 xg cos2 6' 



(33) 



The angle 6 is connected, in both cases, with the collision 
energy in laboratory reference frame (in- vacuum kinemat- 
ics): 



eiab 



2mN 



{kp^ + kp2 — 2A;pi kp2 cos 6). 



(34) 



The exact thermal conductivity of neutrons Kcxact is re- 
lated to the variational thermal conductivity Hvar through 
the correction factor C: 

12 

Koxact = Kvar C'(Ai^), C(Ak) = — iJ(A/f), (35) 

5 



where 



H{Xk) 



3- A 



K 



V 

^ (n + 1 



4n + 5 



(n + l)(2n + 3)[(n + l)(2n + 3) - A 



(36) 



is the function introduced by Sykes & Brooker (1970). 
For the range of particle Fermi wave-numbers fcpn and 
fcpp considered in the present paper, the parameter Xk 
determined from Eq. (HH) varies from 0.79 to 0.97. 
Corresponding values of C{Xk) lie between 1.20 and 1.22. 
It seems reasonable, therefore, to adopt the constant cor- 
rection factor C = 1.2 to the values of Kvar- 
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2.8. Comparison with other works 

Let us compare our results with the results of several pre- 
vious works at lower densities. First of all, we mention two 
papers by Wambach, Ainsworth & Pines ( 1993[ ) (hereafter 
WAP) and by Sedrakian et al. ( |1994[ ). WAP considered 
pure neutron matter at densities about no = 0.16 fm~'^ 
using Landau theory. They calculated momentum depen- 
dent Landau parameters, quasiparticle transition ampli- 
tudes and transport coefficients (thermal conductivity, vis- 
cosity, and spin diffusion coefficient). As follows from their 
Fig. 10, panel 4 (where the quantity kT/c should be ex- 
pressed in MeV/fm^ rather than in K/fm^ as printed), the 
thermal conductivity is k^^^ w 2 x 10^^ ergs s^^ cm^^ 
K"-'^, if we take n„ = and T — 2 MeV. (This particu- 
lar temperature is chosen for comparison with the results 
of Sedrakian et al. ( 1994 ); other authors neglected tem- 
perature effect on the transition probabilities and, conse- 
quently, found the standard dependence of the ther- 
mal conductivity.) The value of the Landau parameter 
at zero momentum transfer from WAP's Fig. 7 yields the 
effective mass m* « 0.9mN at n„ = no. Using our in- 
vacuum formula with the same to* and omitting np colli- 
sions, we obtain k « 2.9 x 10^" ergs s~^ cm~^ K~^. With 
the medium effects included by using the results of Sect. 
2.6 we get about 4.7 x 10^" ergs s^^ cm^^ K^^. In this 
case, we have taken into account rather weak dependence 
of the in- medium correction i^„, Eq. (^0|), on fcpp- 



The calculation by Sedrakian et al. (1994), based on a 
thermodynamic T-matrix approach, yields « 2.6 x 10^*^ 
ergs s""'^ cm"-*^ which agrees surprisingly well with 

our in- vacuum result for the same T=2 MeV and n„ — 
uq. However, the treatment of Sedrakian et al. (1994) is 
limited to rather high temperatures T > 1 MeV at p ~ pq. 
At lower temperatures, the authors observe a divergence 
of their nn scattering cross-sections indicating the onset 
of neutron superfluidit y. It is worth noting also, that Fig. 
2b of Sedrakian et al. (1994) reproduces inaccurately the 
thermal conductivity obtained by WAP. For instance, it 
gives 2.7 X 10^" ergs s~^ cm~^ K~^, for T and n„ under 
discussion, that is 7.5 times smaller than the actual value 
quoted above. 

Finally, our thermal conductivity based on in-vacuum 
cross-sections can be directly compared with the neutron 
thermal conductivity obtained by Flowers & Itoh ( 19791 ). 
Their results are valid for a model of dense asymmetric nu- 
clear matter proposed by Baym et al. (|1971). From Fig. 



1 of Flowers & Itoh (1979) we have < 



10 



23 



ergs s 



2.7 X 10" g cm-3 and T = 10^ K. 



cm ^ K ^ at p 
Corresponding particle number densities are (Baym et al., 
1971| ) n„ w 0.153 fm"^, and Up w 0.006 fm-^. The nucleon 
effective masses were taken equal to their bare masses. In 
this case, our formulae yield the value k„ « 1.8 x 10^^ 
^, which is more than 5 times smaller. 



ergs s ^ cm ^ K 



Note that Flowers & Itoh ( 1979| ) reported good agreement 
between their results obtained in two ways: from the mea- 
sured vacuum phase shifts and in the spirit of Landau 
theory. A qualitative agreement of their results with the 



Landau-theory results by WAP is clear from the above 
discussion. However the accuracy of their vacuum phase 
shift approach is questionable since k^^ disagrees notice- 
ably with our in-vacuum thermal conductivity. The latter 
is virtually exact (see Sect. 2.(;) and should be the same 
as obtained from the measured vacuum phase shifts. 

Summarizing we stress a large scatter of results ob- 
tained using various approaches. In this respect, a reliable 
calculation of the Landau parameters and momentum de- 
pendent quasiparticle amplitudes in the density range (2 
- 3) po, at least for pure neutron and symmetrical nu- 
clear matter, is highly desirable. It is especially true in 
view of relatively large discrepancies at low densities be- 
tween the calculations which employ the in-vacuum cross- 
sections (like ours) and those performed in the frames of 
the Landau theory. 

3. Effects of Superfluidity of Nucleons 

The equations of Sect. |^ give the neutron thermal con- 
ductivity for non-superfluid nucleons. Now we focus on 
the effects of nucleon superfluidity. 



It is well known (e.g., Tilley & Tilley, |1990D that the 
heat transport problem in superfluid matter is compli- 
cated by the appearance of convective counterflow of nor- 
mal part of fluid component in the presence of a temper- 
ature gradient. Macroscopic counterflow can carry heat 
producing an effective thermal conductivity. A study of 
this effective conductivity is a delicate task which is out- 
side the scope of the present work. We will consider the 
effect of superfluidity on the microscopic (diffusive) ther- 
mal conductivity. This problem is also complicated and 
we will adopt some model assumptions. Our considera- 
tion is basically close to that used in the studies of trans- 
port properties in superfluid "^He by Bhattacharyya et al. 



( |1977D and Pethick et al. (|1977|). 

It is gener ally beheved (see, e.g. 



Takatsuka 



Tamagaki, 1993) that the nucleon superfluidity in the neu- 
tron star cores is of the BCS type. The proton pairing oc- 
curs mainly in the singlet ^Sg state. The neutron pairing 
can be of the singlet type only at not too high densities, 
p ^ Po, that is near the outer boundary of the neutron 
star core. At higher densities, the nn interaction in the 
^Sq state becomes repulsive, and the singlet-state neu- 
tron pairing disappears. However, the nn interaction in 
the triplet •^P2 state is attractive, and the neutron super- 
fluidity at densities p ^ po is most likely to occur in this 
state. We consider the singlet-state proton superfluidity, 
and either the singlet-state or triplet-state neutron su- 
perfluidity. 

Microscopically, the superfluidity leads to the appear- 
ance of an energy gap S in the superfluid quasiparticle 
dispersion relation near the Fermi surface (j/c — fepj <C kp): 



sg 



n(OV^^ 



(37) 



where ^ = hvp{k—kp) is the normal-state quasiparticle en- 
ergy measured with respect to the chemical potential and 
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v-p is the normal-state quasiparticle velocity at the Fermi 
surface. In the case of the singlet-state pairing, the energy 
gap (5 = (5o is isotropic, i.e., independent of orientation of 
the particle momenta with respect to a quantization axis. 
The temperature dependence of Sq can be fitted as (e.g. 
Levenfish & Yakovlev, 



1994b) 



y 



T 



^71^(1.456 -^ + 1:1^), (38) 



where r = T/Tc and Tc is the critical temperature. 

In contrast to the isotropic singlet-state pairing, the 
pairing in a triplet-state produces an anisotropic gap 
which depends on orientation of nucleon momentum. If 
the direction of the quantization axis were fixed for some 
physical reason, this would complicate the analysis be- 
cause the conductivity would be anisotropic with respect 
to the quantization axis. For the sake of simplicity, it 
is usually assumed that the orientation of the quantiza- 
tion axis in the neutron star cores is uncorrelated with 
rotational axis, temperature gradient, or magnetic field. 
Therefore, matter of the neutron star cores is treated as a 
collection of microscopic domains with arbitrary orienta- 
tions of the quantization axis. 

Thus we assume that the diffusive thermal conductiv- 
ity is isotropic even in the presence of the triplet-state 
neutron pairing. Under this assumption, we may use an 
effective dispersion relation with an isotropic energy gap 
S = Si which greatly simplifies our analysis. It would be in- 
appropriate, though, to employ the same temperature de- 
pendence of the gap as for the singlet-state pairing gap Sq, 
Eq. (38). For instance, this would yield inaccurate gap pa- 
rameters at T ^ Tc- We will take advantage of the detailed 
analysis by Levenfish & Yakovlev (1994a), Levenfish & 
Yakovl ev (|1994bD , Yakovlev & Levenfish ( |1995| ), Yakovlev 
et al. ( |1999| ) of the nucleon heat capacity and various 
neutrino-energy loss rates in the neutron star cores for 
the case of triplet-state pairing of the neutrons with zero 
projection of the total angular momentum of Cooper pairs 
onto the quantization axis (mj — 0). According to these 
studies, one can obtain qualitatively accurate results us- 
ing the isotropic gap Si equal to the minimum value of the 
angle-dependent gap on the Fermi surface. The temper- 
ature dependence of this artificially i sotro pic gap can be 
fitted as (e.g., Yakovlev & Levenfish, 

Si(T) n / 1.188' 



19951 ) 



T 



0.7893 



(39) 



Low-lying excited states of a nucleon superfluid corre- 
spond to the presence of a dilute gas of superfiuid quasi- 
particles. The time and space evolution of the quasipar- 
ticle distribution function is governed by the transport 
equation. On the left-hand-side (streaming terms) there 
appear two terms containing derivatives of the superfluid 
gap but they cancel out. As a result, the left-hand-side 
acquires the standard form, identical to that for normal 
quasiparticles. 

However, the collision term is, in general, more compli- 
cated than that for normal nucleon liquid. The scattering 



amplitude for quasiparticles in the superfluid is a linear 
combination of the normal-state amplitudes. These effects 
were considered in detail by Bhattacharyya et al. (1977) 
at T ^ Tc and by Pethick et al. ( |1977D at T < Tc for liquid 
■^He. The accurate expressions for the superfluid scatter- 
ing amplitudes are complicated. In view of the large uncer- 
tainty in the amplitudes even in the case of non-superfluid 
matter, we will keep the transition probabilities the same 
as in the normal case. Thus we restrict ourselves to the 
most important phase space effects stemming from the 
modification of the quasiparticle dispersion relations and 
resulting in a specific temperature dependence of the col- 
lision frequencies. Our results at T ^ Tc and T ^ Tc are 
in q ualita tive agreement with those of Bhattacharyya et 
al. ( |1977| ) and Pethick et al. (|l977| ). 

In order to calculate the diffusive thermal conductiv- 
ity in superfiuid matter we adopt the same variational ap- 
proach which was used in Sect. ^ but with the modified dis- 
persion relations, Eq. (37), containing energy gaps. In par- 
ticular, the deviation function is given by the same Eq. (^, 
where the velocity is now given by u = i'f|CI/vs+^- 
Thus, the velocity of superfluid particles varies sharply 
near the Fermi surface and vanishes ai k — kp. 

The temperature dependences of the collision integrals 
can be obtained using the expressions for the energy inte- 
grals, Ja, but with dispersion relations of superfluid quasi- 
particles. Notice that the combination of Fermi-Dirac dis- 
tributions in the collision integral remains a sharp function 
of momenta which justifies our energy-angular decompo- 
sition (Sect. 2.3). Then from Eq. (fiSh we obtain 



7 - 7(°)7^ 



(40) 



where Jq°^ is an energy integral ( |l6| ) in non-superfluid 
matter and TZa describes the effect of nucleon superfluid- 
ity. As we use the same Q12 as in the normal state, the fac- 
tor TZa represents, in our approximation, the overall effect 
of nucleon superfluidity on the collision integral. Our anal- 
ysis shows that TZa < 1 in the presence of neutron and/or 
proton superfluidity. In other words, the superfluidity al- 
ways suppresses the nn and np collision frequencies. There 
is no discontinuity of TZa when a superfluidity is switched 
on with decreasing T. The suppression becomes very large 
{TZa < 1) for strong superfluidity ( T > Tc ). 

Following Levenfish & Yakovlev (1994a) and Yakovlev 
& Levenfish ( |1995D we introduce the notations 

^ = J^^ y^J,, z = ^^sgn{x)y/^P~^^. (41) 

Then the correction factors for the np collisions can be 
written as: 



T^pi{yn,yp) = TTl I I I I dxndxpdx',^dx' 



T^p2{.yn,yp) 



27T^J J J d^"dXp..„.u,p 
X S{Zn + Zp + z'^ -f Zp) 
X f{Zn)fiZp)f{z',^)f{Zp)xl 

15 
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X S{Zn + Zp + z'j + z'p) 

X f{Zn) f{Zp) /«) fiz'p) {-Xr,x'J, (42) 

where /(z) = (e^ + l)~^ = l—/(—z). The reduction factors 
for the nn colhsions are immediately expressed through 
7? • 



(43) 



The normalization factors in Eqs. (^2|) are chosen to sat- 
isfy the condition TZa ~ 1 f^or vanishing gaps. Furthermore, 
it is possible to analyze the asymptotic behaviour of TZa 
in the cases in which the neutron and/or proton super- 
fluidity is strong (y 3> 1). This can be done in the stan- 
dard manne r (cf. L evenfish & Yakovlev, 1994a , Levenfish 
& Yak ovlev , |l994b| , Yakovlev & Levenfish, |l995| , Yakovlev 
et al., 199£ ) Let us summarize the results. 

If the neutron superfluidity is strong (j/„ 2> 1), the nn 
collision frequency is affected by the factor 



15 



-2y„ 



(44) 



In addition, we have calculated TZnaiHn) numerically 
and proposed analytic fits which reproduce the numerical 
results and the above asymptotes: 



2 r 1 

- 0.9468 + V(0-0532)2 + 0.5346 2/2 

O 

X exp 



0.377- V(0.377)2 -1-4 y2 



i(l + 1.351 y2)2 



X exp 

1 r 



0.169 - v/(0.169)2-h9y2 



2 

X exp 
1 

+ 2 



0.6242 + \/(0.3758)2 -f 0.07198 y2 



3.6724 - v/(3.6724)2 + 42/2 
(1 + 0.01211 y2)9 



X exp 



7.5351 - V(7-5351)2 -H9y2 



(45) 



Maximum errors of these fits do not exceed 0.5%. 

Regarding np collisions, we have calculated 
TZpaiymUp) numerically for the cases in which ei- 
ther neutron or proton superfluidity is absent and fitted 
the numerical results by the expressions: 



Rpi{yn,0) 

Rp2{yn,0) 

T^pi{o,yp) 



0.4459 + V(0-5541)2 -I- 0.03016 ?;2 
X exp [2.1178- v/(2-1178)2-hy2 



0.801 + V(0-199)2 + 0.04645 y2 



X exp 
1 



2 

X exp 



2.3569- V(2.3569)2 
0.3695- 
2.4451 



(0.6305)2 0.01064 y2 



(2.4451)2 y2 



(l + 0.1917y2)i-4 



Tip2{Q,yp) 



X exp 
0.0436 



4.6627- W (4.6627)2 + 42/2 



(4.345)2 -^19.55 2/2-3.345 



X exp 
+ 0.0654 exp 
+ 0.891 exp 



2.0247- ^/(2.0247)2 -K2/2 



8.992 



(8.992)2 -f 1.5 2/2 



9.627- w'(9.627)2 + 92/2 .(46) 



These fits reproduce numerical results calculated in the 
range ?/„ < 30 or 2/p < 30. The maximum fit error does 
not exceed 1%, 2%, 0.3%, 1.4%, respectively, and the mean 
fit error is about twice lower. 

The most difficult is the case of np collisions in matter 
in which neutrons and protons are superfluid at once. In 
this case we have calculated 7?.pi. 2(2/71, 2/p) numerically on 
a dense grid of 2/n and 2/p (2/n < 12, 2/p ^ 12) and obtained 
analytic fits to these results: 



Tlpi{yn,yp) 



(0.7751 + 0.4823 u„ + 0.1124 Up 



0.04991 ul + 0.08513 u^Up + 0.01284 ulup) 
g-2u+ (0.2249 -f 0.3539 u+- 0.2189 M_ 
0.6069 u„u_ -I- 0.7362 UpU+) , 



(47) 



where up = ^yj^ + (1.485)2 - 1.485, (3 = -h, -, n, p, 

2/_ = min(2/„,2/p) and y+ = max(2/„, 2/p). The mean fit 
error is about 3%, and the maximum error of 14% takes 
place at 2/n ~ 10 and 2/p ~ 12, where TZpi is exponentially 
small. 

The fit to Tlp2 reads 

'R-p2{yn, yp) = e-"+-"- (1.1032 -I- 0.8645 m„ + 0.2042 Up 
+ 0.07937 + 0.1451 w„Wp + 0.01333 ulup) 



-2tt4 



-0.1032 - 0.2340 u+ -f 0.06152 u„7i4 



0.7533 M„u_ — 1.007 MpM4.) 



(48) 



where up = ^y^ -H (i.Ytjij' — i.YDi witn p = +, —, n, p. 

The mean fit error is about 3.6%, and the maximum error 
of 14% again takes place at 2/n ~ 10 and 2/p ~ 12 where 
7?,p2 is exponentially small. The accuracy of our fit for- 
mulae, Eqs. ( ^ ) and (^8|), is quite sufficient for practical 
applications. One can easily see that these fit formulae 
are in good agreement with Eqs. (^5|) at yn — yp and with 
Eqs. ( ^6| ) either at 2/n = or at 2/p = 0. 

The neutron superfluidity affects also the left-hand 
side of the kinetic equation given by the integral (|l^) . For 
a non-superfluid matter, we obtained C„o = k^^T/{im'^) 



(Sect. 2.2). In the presence of the neutron superfluidity, 
we can write 



Cn{yn) = CnoTZciyn), 



^c(2/n) 



6 



(e-. - 1)2 



(49) 
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In addition, the neutron superfluidity affects the neutron 
heat flux in Eq. (|l|). It turns out that this effect is de- 
scribed by the same factor TZc- Consequently, instead of 
Eqs. ^ and (|l2|) we can express the diffusive thermal 
conductivity in the form 



3m,* 

Sn2 Kn2 'R-n2{yn) 



7^r 



'12, 



(50) 



p2 Tip2{yn,yp) 



^p2 



1 



2 ^pi Kpi [37^pl(y„ 



,yp) - 'R-p2{yn,yp)] 



The new factor TZciVn) is analyzed easily in the stan- 
dard manner. It is equal to 1 in the absence of superflu- 
idity, and behaves as TZciyn) ~ 3^2 (j/„/7r)'^/^ exp(— ?;„) 
in the limit j/„ ^ oo. We have calculated this factor nu- 
merically on a dense grid of y„. These results and the 
asymptotes are accurately fitted as 



7^c(yn) = 0.647+ V(0.353)2 + 0.109 y2 



1.5 



X exp 1.39- V(l-39)2 



The maximum fit error does not exceed 1%. 



(51) 



4. Results and Discussion 

Let us summarize the results of Sects. ^ and ^and present 
practical formulae for the neutron thermal conductivity /t„ 
in the neutron star cores: 

_ 7r2fc|rn„T„7^c(2/n)C 



3m* 

7.2 X 1023 Ts7^^(y, 
IQi^s-i 



m* 



no 



ergs cm ^ s "'^ K ^ . 



(52) 



Here, uq — 0.16 fm"^ is the normal nucleon density, Tg 
is temperature in units of lO'* K, is the Boltzmann 
constant (now presented explicitly), C=1.2 is the correc- 
tion factor discussed in Sect. 2/7; TZciVn) is the superfluid 
reduction factor. The latter factor depends on the neu- 
tron gap parameter yn and is given by Eq. (^l]) . The gap 
parameter ?/„ is determined either by Eq. (|38|) for singlet- 
state neutron pairing or by Eq. (p£ ) for triplet-state neu- 
tron pairing as explained in Sect. |. Let us remind that 
we consider the diffusive thermal conductivity and do not 
analyse convective heat transport in superfluid matter. 

According to Eqs. ( ^l]) and (^^, the nn and np colli- 
sion frequencies in Eq. (|52|) are 



3.48 X 10 



15 



m^ 

nir, 



[^n2 ^n2 '^«2(2/«) 

3^1"^ Kni [7^„l(2/„) - 7^„2(y„)]} s-\ 



26m;m;2fc2T2 



3.48 X 10 



15 



m^ 



^ Ti 



{^^^2 Kp2 'R-p2{yn,yp) +0.5Kpi 

S^^i [37^pl(2/„,yp) - 7^p2(yn,2;p)]} 



(53) 



Here, S^i"' ( a = nl, n2, pi, p2) is a normalized transition 
probability integrated over momentum transfers of collid- 
ing particles (Sects. 2.4-p]6) for free nucleons [expressed in 
mb and fitted by Eqs. describes the in-medium 

effects on the squared matrix element [the fits are given 
by Eqs. (^0|)]; TZa is a corresponding superfluid reduction 
factor whose fits are given by Eqs. (^)-(^8|). Thus we ob- 
tained a simple description of the diffusive neutron con- 
ductivity in the density range from 0.5/3o to 3/9o (see Sect. 
H) . Although we have employed the Bonn model of AW in- 
teraction, we hope that the results could be used, at least 
semi-quantitatively, for a wide class of equations of state 
if one takes effective nucleon masses m^, Fermi momenta 
fcpAf, and superfluid gaps appropriate to these equations 
of state. Moreover, our results are presented in the form 
which enables one to include the in-medium effects from 
more elaborated (future) calculations of the nucleon quasi- 
particle transition probabilities. It would be sufficient to 
recalculate four in-medium correction factors Ka [defined 
in Eq. ^] while aU other quantities in Eqs. (|2|)-(|3|) 
will remain untouched. 

The results are illustrated in Figs. 0-^. Figure |^ 
presents the neutron thermal conductivity versus density 
at T = 5 X 10^ K for non-superfluid matter. Here, we 
use the simplified model of free degenerate ne utron s, pro- 
tons and electrons (e.g., Shapiro & Teukolsky, 1983). This 
model yields an extremely soft equation of state, with very 
low fraction (about 0.6% a,i p — po) of protons and elec- 
trons. The effective masses of nucleons are set equal to 
their bare masses, and no in-medium effects are included 
in the nn and np collision frequencies. The solid curve de- 
picts the full neutron conductivity, which takes account of 
nn and np collisions. The dashed line shows the contribu- 
tion from nn collisions alone. The nn and np collisions are 
seen to be of comparable importance despite the extreme 
smallness of the proton fraction. This is related to the fact 
that the np collisions occur at smaller laboratory energies 
and scattering angles than the nn ones, while the NN dif- 
ferential cross sections tend to grow with decrease of both 
these parameters. We present also the electron conductiv- 
ity (the dott ed line ) calculated using the results of Gnedin 
& Yakovlev (|l995| ). 

Figure |^ demonstrates the significance of in-medium 
effects. It displays the ratios of in-medium to in-vacuum 
nN collision frequencies versus Fermi wavenumber of neu- 
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Fig. 1. Thermal conductivity of neutrons and electrons 
versus mass density (in units of po = 2.8 x 10^'' g cm^'^) 
of non-superfluid npe matter in the model of free particles 
at T = 5 X 10^ K; = tun. 




Fig. 2. Ratios of in- medium to in- vacuum nn and np colli- 
sion frequencies versus fcpn at A:fp= 0.3 fm~^ (solid lines), 
0.6 fm^i (dashes) and 0.9 fm^^ (dots); = 0.8 m,N- 



trons fcpn (notice that fcpn ~ 1-7 fm~-^ at p — p^) at three 
values of the proton Fermi wavenumber, k-p-p = 0.3, 0.6 
and 0.9 fm~^. Let us emphasize again that we use one 
particular model of AW interaction in medium which is 
valid, strictly speaking, for symmetric nuclear matter as 



Fig. 3. Thermal conductivity of neutrons and electrons 
versus mass density in non-superfluid npe matter with a 
moderately stiff equation of state (Prakash et al., 198S 
Page & Applegate, |1992D, = O.Tto^v at T = 5 x 10 



K. Solid curve is calculated with the in- medium nn and np 
scattering cross sections while other neutron-conductivity 
curves are obtained with the in-vacuum cross sections. 



described in Sect. |2.6| . The in- medium effects are seen to 
be significant, especially for np collisions, and determined 
mostly by fcpn. For instance, at fcp„ « 2 fm^^ the in- 
medium effects reduce Vnp by a factor of 3-6. The depen- 
dence of the ratios in question on the proton fraction (on 
fcpp) appears rather weak. 

Figure |^ shows the neutron and electron thermal con- 
ductivities versus density at T = 5 x 10^ K for non- 
superfluid matter described by a moderately stiff equa- 
tion of state proposed by Prakash et al. ( 1988 ) [the ver- 
sion with the compression modulus Ka — 180 MeV, and 
t he sa me model function F{u) as in Page & Applegate 
(1992)]. This equation of state is much more realistic than 
the model of free particles employed in Fig. [l| it yields 
much larger fractions of protons and electrons (e.g., about 
6% ai p = po). The effective nucleon masses are put equal 
to 0.7 mj^. The solid curve takes into account the con- 
tribution from nn and np collision frequencies evaluated 
with the in- medium cross sections (Fig. ^). The curve is 
plotted up to yO « 3 po; in the density range where the in- 
medium scattering cross sections are available (Sect. ^.61 ). 
Other neutron-conductivity curves are calculated with the 
in-vacuum cross sections [Ka — 1 in Eqs. (|5^)]. In con- 
trast to Fig. 1^, the electron thermal conductivity exceeds 
the neutron one at all densities, which is explained by the 
much higher fraction of electrons for the adopted equation 
of state. Comparing Figs. |l| and ^one can deduce that the 
electron conductivity is more sensitive to the equation of 
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Fig. 4. Diffusive thermal conductivity of neutrons (solid 
lines) and electrons (dashes) versus temperature at p = 
2po for the same equation of state as in Fig. |[ The curves 
marked SFO refer to the non-superfluid case; the curves 
marked SFl imply = 3 x 10*^ K, T^p = 10^ K, while 
the curve SF2 corresponds to Tc„ = 4 x 10^ K, T^p = 10^ 
K (singlet-state proton and triplet-state neutron super- 
fluidities). The in- medium effects in the scattering cross 
sections are not included. The dashed curve SFl is equally 
valid for the case SF2. 

state than the neutron one. This is because the number 
density of electrons varies more strongly with the equation 
of state than the number density of neutrons. In fact, the 
variation of the neutron conductivity is mostly caused by 
the difference in effective nucleon masses between the two 
models. 

Finally, Fig. || demonstrates the effects of neutron and 
proton superfluidity. It shows the temperature dependence 
of the diffusive thermal conductivity for the same model 
of matter as in Fig. ^ at p = 2po. Solid and dashed lines 
display the neutron and electron thermal conductivity, re- 
spectively. The curves marked as SFO correspond to non- 
superfluid matter; the curves SFl correspond to the model 
in which = 3 x 10^ K, T^p = 10^ K , i.e., the proton 
superfluidity is stronger than the neutron one. Finally, the 



curve SF2 is for the case 



4 X lO'^ K, T„ 



10" K, 



in which the neutron superfluidity is stronger. The pro- 
ton superfluidity is assumed to be of singlet type [the de- 
pendence of the gap on the temperature is given by Eq. 
(^)] while the superfluidity of neutrons is taken to be of 
triplet type with the gap made artificially independent of 
the direction of particle momentum [the temperature de- 
pendence of the gap is given by Eq. (|3^)]. Note that the 
electron thermal conductivity is independent of superfluid 
state of neutrons and is the same for models SFl and SF2. 



All the solid curves are obtained with in- vacuum cross sec- 
tions. In the absence of superfluidity, both conductivities, 
K„ and Ke, behave as T~^. When temperature falls be- 
low Ten or Tcp, the conductivity depends on T in a more 
complicated manner. 

If Ten > Tcp (curve SF2, dashed curve SFl is equally 
valid for the superfluid model SF2), and the tempera- 
ture falls below Ten the electron conductivity remains un- 
affected but the diffusive neutron thermal conductivity 
drops exponentially due to sharp decrease of the density of 
heat carriers — neutron quasiparticles whose momenta are 
sufficiently close to the neutron Fermi momentum. When 
T falls below Tcp the electron thermal conductivity under- 
goes a slight reduction associated mainly with variation in 
the plasma screening length (the proton screening which 
is efficient at T > Tcp dies out leaving alone the weaker 
electron screening). 

If Tcp > Ten (curves SFl), the proton superfluidity oc- 
curs first with decreasing T. It switches off the np collisions 
and thus enhances k„. When the temperature reaches Ten, 
the neutron conductivity continues to grow over the non- 
superfiuid one. At further decrease of T the neutron con- 
ductivity is determined by competition of two factors: the 
suppression due to the decline of the density of the heat 
carriers, and the enhancement due to lowering the nn colli- 
sion frequency (while the np collisions are almost switched 
off). If T ^ Ten both factors nearly cancel each other out, 
and the conductivity, although deviates somewhat from 
the non-superfluid one, but reproduces the standard tem- 
perature dependence T^^. Analogous effect was demon- 
strated by Pethick et al. (1977) in the case of superfluid 
^He at T < Tc. 

We see that the nucleon superfluidity affects the neu- 
tron thermal conductivity in different ways, producing ei- 
ther enhancement or suppression, which is incorporated 
easily in computations using Eqs. ( |5^ and (53). 

We have considered thermal conductivity of matter 
composed of neutrons, protons, electrons, and muons at 
p Pa. Let us discuss briefly thermal conductivity in 
hyperonic matter which may appear at higher p (e.g., 
Balberg et al. 1999). It may contain S~, A and other hy- 
perons. Theoretical knowledge of hyperon interaction is 
rather uncertain but numerous models of hyperonic mat- 
ter give qualitatively the same results. S" and A hyperons 
appear at about the same density. Their relative fractions 
increase sharply with growing p while the fractions of elec- 
trons and muons become much lower than in the absence 
of hyperons. The fraction of A hyperons is typically larger 
than that of E~ hyperons. Qualitatively, electrons and 
muons are replaced by hyperons while neutrons are 
are partially replaced by A hypeons. Both replacements 
have important consequences for thermal conduction. 

It is likely that thermal conductivities limited by 
Coulomb and nuclear interactions will remain decoupled. 
Lowering the electron and muon fractions reduces the ef- 
ficiency of the electron and muon Coulomb thermal con- 
ductivities. These conductivities can be de termin ed easily 
using the results of Gnedin & Yakovlev (1995). On the 
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other hand, the efficiency of thermal conductivity, limited 
by strong interactions, becomes higher. Neutrons and A 
hyperons may become the main heat carriers of compa- 
rable importance. Their thermal conductivity will be de- 
termined by collisions of these particles with themselves 
and with other strongly interacting particles (protons, , 
etc.). Their effective collision frequencies are likely to be 
coupled. The theory of thermal conductivity of neutrons 
and A hyperons may be constructed using the formalism 
of Flowers & Itoh (|1979|) and Anderson et al. ([1987D. The 



main problem would be to calculate the scattering cross 
sections of strongly interacting particles of various species 
taking into account many-body effects. Hyperons, like nu- 
cleons, may be in superfluid state (e.g., Balberg et al. 



199£) forming multicomponent superfluid. The effects of 



superfluidity on the diffusive thermal conductivity can be 
included in the same manner as in Sect. 3. Consideration 
of convective heat transport would be even more com- 
plicated than in nucleon matter (see above). Therefore, 
thermal conduction in hypcronic matter remains an open 
problem. It deserves a separate study which goes far be- 
yond the scope of present paper. 

5. Conclusions 

We have calculated the thermal conductivity of neutrons 
due to nn and np collisions in the neutron star cores. The 
results are valid for densities p from O.Spo to 3po (see Sect. 
^). We have included the effects of possible superfluidity 
of neutrons and protons (Sect. ^) and obtained practical 
expressions (Sect. ^) for the diffusive thermal conductiv- 
ity valid for a wide class of models of dense matter. The 
results can be generalized to the case in which hyperons 
are present in dense matter along with nucleons. 

Our results, combined with those of Gnedin & 



Yakovlev ( 1995 ) for the electron thermal conductivity, 
provide a description of the transport properties in the 
neutron star cores. This thermal conductivity is needed 
for numerical simulations of cooling of young neutron stars 
with non-isothermal core (age up to 10 — 100 years). 
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